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Flavor dependence of normalization constant for an infrared renormalon 
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An ansatz is proposed for the flavor dependence of the normalization constant for the first IR 
renormalon in heavy quark pole mass. 


The ultraviolet (UV) and infrared (IR) renormalons in field theory are fascinating objects that give rise to factorially 
growing large order behavior in perturbation theory. In quantum chromodynamics (QCD) the IR renormalons can 
be understood within the framework of operator product expansion (OPE) [H-Q. An IR renormalon in a Wilson 
coefficient in the OPE causes ambiguity in the resummed asymptotic series of the Wilson coefficient, which is to 
be cancelled by the vacuum condensate of an higher dimensional operator. This cancellation of the ambiguities in 
the resummation and vacuum condensate determines the nature of the renormalon singularities in the Borel plane. 
Specifically, the functional form of the condensate in the strong coupling Us, which can be determined up to an 
overall constant by renormalization group (RG) equation, determines the power of singularity as well as its location 
0. However, the residue of the singularity, which is the normalization constant of the renormalon-caused large order 
behavior, is not known, but it can be calculated by a perturbation method @, 0. While the normalization can be 
expressed in a convergent series, with finite order perturbation it can only be calculated approximately. Its exact 
form is a nonperturbative quantity and so far there is no known way to calculate it. It would thus be very interesting 
if we somehow find the exact form for the normalization. 

Our purpose in this Letter is to present an ansatz for the normalization for the first IR renormalon in the heavy 
quark pole mass. In this case the perturbation method for the normalization yields a series that converges rather 
quickly, which allows to determine the normalization within a few percent of uncertainty using the known first three 
perturbative coefficients ©[HI- The reason for this rapid convergence and accuracy may lie with following two facts. 
First, the renormalon singularity is relatively soft, with its singularity of (I — 26)“G+/3i/2/3o)^ fQj- example, compared 
to that of the first IR renormalon in Adler function, which is of (I — where /3o,/3i are the first two 

coefficients of the beta function and b is the complex variable for the Borel plane. With vanishing flavor number 
Nf = 0, for example, the singularities are (1 — and (I — 6/2)“^ ®®, respectively. Secondly, the IR renormalon 

is the closest singularity to the origin in the Borel plane. This condition is important because the normalization in the 
perturbation method is evaluated on the boundary of convergence disk and for that to work the singularity should be 
closer to the origin than any other singularities. Of course, any renormalon singularity can be moved by a conformal 
mapping to be the closest one to the origin, this step is not required with the first IR renormalon for the pole mass, 
because it is already the closest one, and this seems to help the convergence. The normalization constant obtained 
from perturbation method was confirmed by recent lattice calculation of large order behaviour of static energy 

To get the ansatz let us assume that the condensate in the OPE is of dimension n and is proportional to 

A5cD-^e-^a”^(l + 0(a.)) (1) 

and the associated renormalon singularity is of the form: 


(I - 0/00^+- ■ ' ' 

Then the imaginary part, which is ambiguous, of the Borel integral 

= ±JV sm{i'7T)r{-iy){bo/a~'' {1 + 0(q!s)) (3) 

is to be cancelled by ambiguity of the form (HD. Thus identifying ([3|) with CD we get 

(jr 2/3o \v 

bo = nl2, zz = n/3i/2/3g , J\f = — — ^ -r, (4) 

sin(z/7r)l {—ly) 

where C is an unknown proportionality constant. Now putting 

C = Co/(0 (5) 
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Nf 

0 

1 

2 

3 

4 

5 

A/pert. 

0.6081 

0.5981 

0.5867 

0.5730 

0.5551 

0.5288 

W/Wpert. 1.0047 

0.9979 

0.9973 

1.0081 

1.0403 

1.116 


TABLE I. 


Comparison of the ansatz with normalization from perturbation method for varying number of flavors. 


with /(O) = 1, where Co is determined so that the normalization agrees with that from the large-/3o approximation 
at jz = 0, we get 


, cofHC-^r 

sin(jz7r)r(—zz) 


( 6 ) 


We shall now consider the normalization for the pole mass, for which n = 1. While in this case the ambiguity in 
Borel summation is not canceled by a condensate of dimension one operator but by the static inter-quark potential 
M, the ambiguity is nevertheless of the form (|^ and the relations (|1]) are still valid. For this case the large-/Sg 
approximation fixes Cq in the MS scheme as |1^17l | 


Co = , 


( 7 ) 


where Cp = 4/3 for QCD. Of course f{v) is not known, and we can only guess it by comparing the normalization 
with values from the perturbation method. This gives 




( 8 ) 


and, finally, 


Cpeii^r 

sin(zz7r)r(—iz) ’ 


( 9 ) 


for the MS scheme, where zz = /3 i/2/?q. This is our proposed ansatz. 

How well does this ansatz work? Table |I] compares the ansatz with the normalization from the perturbative method 
for varying number of flavors, the details of which will be discussed shortly. The agreement is quite impressive for 
up to iVy = 4. It is to be noted that the perturbative method tend to work better for smaller Nf. Of course, this 
does not mean that the ansatz is necessarily correct. One test may be to expand ([9]) for small zz and compare it with 
subleading corrections in large-/3o approximation, for which only partial result exists [l8l| . Expanding at zz = 0, 
we have 

A/" = —[1 -b (te -b log(/3o/2))zz -b C(zz2)] ^ 

TT 

The Euler constant term agrees with the corresponding term in [l8l| . It is worth noting that the two transcendental 
numbers tt and 7 e in OT both arise from 

sin(zz7r)r(—zz) (11) 


in the denominator in ([2]). This is a strong indication that the term should be, at least, part of the exact form for the 
normalization. 

It may be tempting to use the same idea to make an ansatz for other renormalons, for example, such as the first IR 
renormalon in the Adler function, which is associated with the gluon condensate, for which n = 4. However, there is 
no accurate estimate of the normalization for the Adler function; The convergence from the perturbative method is 
not fast M- Without accurate numerical estimates for the normalization it is impossible to infer /(zz) for the Adler 
function and that blocks such an attempt. Nevertheless, if we assume (|S]) is applicable to the Adler function as well, 
then we have an ansatz for the hrst IR renormalon of the Adler function 


Af=- 


^Cpefifr 


sin(zz7r)r(—zz) ’ 

where zz = 2/3i//3o, and the following large-/3o result is used [13 Haim: 


Co = ^Cpe- 


( 12 ) 


( 13 ) 
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It is interesting to see how the ansatz compares with perturbative results, even though the convergence for them is 
not that good and uncertainty is too large for accurate comparison. From [l^, which uses five-loop Adler function, 
the normalization from perturbation is given as 0.287,0.251,0.208,0.154 for Nf = 0,1,2,3, respectively, and the 
corresponding numbers from the ansatz are 0.327, 0.300,0.282,0.277. As expected the ratios of corresponding numbers 
are not close to unity, but what is remarkable is that numbers from the ansatz and perturbation fall in the same 
ballpark. This is not so a trivial point because there is an order of magnitude difference between the large-/3o result 

— = 8.9 (14) 

TT 

and the numbers from the perturbative method, and there is no obvious reason that a function of such complex form 
as m should give values that are of same magnitude as the perturbative numbers. 

We now show how A/'pert. in Table U were obtained. The bilocal expansion, which interpolates the two expansions 
about the origin and the renormalon singularity, of the Borel transform ?7i(6) for the pole mass is given in the form 

Pil: 




with which the Borel summed pole mass tubr is given by 


TOBR = TO¬ 


MS 


1 + Re(^-^J^ 


where as = as (to^) . The perturbative form of the Borel transform is 


n=0 


which gives the perturbative expansion of the pole mass: 

^pole = ^ ^ PnC^s ) 5 

n—0 

where the hrst three coefficients are given as (23l - l^ 

po = 0.4244, Pi = 1.3621 - 0.1055Af/, p2 = 6.1404 - 0.8597Af/ -k OmilN] . 
The normalization M is given by 

V=i?(l/2), 

where 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 


R{b) = m{h){l - (21) 

and expanding R{b) at the origin using (flTl) the normalization can be evaluated perturbatively [1,0. While this yields 
a convergent series for the normalization it does not fully exploit the bilocal expansion (1151) . especially the expansion 
about the renormalon singularity, of which the coefficients Ci are entirely dependent on the beta function |18l |. and 
with four-loop beta function the hrst two coefficients ci,C 2 are known |10l |. To utilize this expansion we write R{b) in 
a truncated form: 

2 

(l-26)i+"+V[l + ^c,(l-26)*], (22) 

and determine hg, hi, A/" by demanding (|2^ and (1^ with (flTl) give identical expansion about the origin to 0{b‘^). To 
get the numbers in Table U this procedure was performed not in the 5—plane but in the conformally mapped z-plane 
dehned by 

b 

TTh’ 



z = 


(23) 
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where the main advantage of this mapping is that the first UV renormalon at 6 = —1 is mapped away to infinity. The 
first IR renormalon is now at z = 1/3 and the bilocal expansion of R in z-plane is given by 


R{b{z)) = 


E- 

.n=0 


(l_3^)i+-+^[l + ^cf(l-3zr] 


(24) 


i=l 


from which jVpert. in Table |T] were obtained. Note that cf, C 2 are given in terms of ci, C 2 by 


Cl = i^ci , C2 = -C2 - -Cl . 


(25) 


In summary, we have proposed an ansatz for the normalization constant for the first IR renormalon in the pole 
mass and presented an argument that, at least, (HU should be an integral part of the normalization constant. 
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